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SUMMARY 


Silicon pn junction solar cells made with low- 
resistivity substrates show poorer performance than 
traditional theory predicts. The purpose of the research 
sponsored by grant NSG-3018 is to identify and characterize 
the physical mechanisms responsible for this discrepancy. 
Attention concentrates on the open-circuit voltage in 
shallow-j unction cells of 0.1 ohm-cm substrate resistivity. 

For such cells we have considered from a theoretical 
vantage point the multitude of possible mechanisms that can 
occur in silicon devices. We have found that two mechanisms 
are likely to be of main importance in explaining the 
observed low values of open-circuit voltage; (1) recombination 
losses associated with defects introduced during junction 
formation; and (2) inhomogeneity of defects and impurities 
across the area of the cell. To explore these theoretical 
anticipations, we have designed and fabricated various 
diode test structures and have constructed measurement con- 
figurations for characterizing the defect properties and 
the areal inhomogeneity. Preliminary experimental findings 
were made during the period pf the grant, and experimental 
studies now continue under a renewal of grant support. 
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I . INTRODUCTION 


In the twenty years since the first realization of a 
silicon solar cell in 1954 by Chapin, Fuller and Pearson 
at the Bell Telephone Laboratories [1] , substantial research, 
development and manufacturing efforts have been spent to im- 
prove efficiency [2,3] and radiation resistance [4] for space 
applications. The design of commercial silicon solar cells 
for space applications has been stabilized to n on p, 10 ohm- 
cm silicon base material of 2x2 cm area. At present, the 

most advanced commercial solar cells now in production yield 

2 

about 13% efficiency, 39- mA/cm short circuit current and 

2 

585 mV open circuit voltage under 140 mW/cm AMO solar illumina- 
tion. 

Shockley and Queisser [51 have shown that the intrinsic 
limits of the silicon solar cell efficiency is about 30% with 
an open-circuit voltage of about 840 mV. In the idealized 
case they proposed, all loss mechanisms are neglected except 
the intrinsic loss due to interband radiation recombination 
of electron-hole pairs. A semi- empirical extrapolation by 
Shockley and Queisser gave an extrinsic efficiency of about 
22%. Taking into account the various losses more quantita- 
tively, other workers have provided different estimates of 
the maximum efficiency: 22% by Wolf [2] and 18% by Brand- 

horst [3] . Both predict a corresponding open-circuit voltage 
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of about 710 mV when the silicon substrate resistivity is 
reduced below 0.1 ohm-cm. 

Hence comparison of these estimates against the per- 
formance of commercial solar cells demonstrates a disagree- 
ment of nearly a factor of two in efficiency and of 100 to 
150 mV in open- circuit voltage. Moreover, these theoretical 
limits have yet to be achieved in any experimental devices. 

For example, recent work shows that the open-circuit voltage 
saturates at about 600 mV as substrate resistivity is reduced 
to 0.1 ohm-cm [3]. 

This study seeks understanding of the physical mechanisms 
responsible for these discrepancies between experimental 
observation and theoretical estimate. In particular, to the 
present, we have concentrated most effort toward understanding 
the open-circuit voltage in shallow- junction cells with 0.1 
ohm-cm substrate resistivity. The interest here derives from 
three considerations: first, the curve fill factor and the 

short-circuit current are thought to show better agreement 
between theory and experiment for these cells than does the 
open-circuit voltage [6,7]; second, it is believed [6] that 
much larger improvements are possible in open- circuit voltage 
than in either fill factor or short-circuit current; and, 
third, experimentation [6,7] suggests that ' performance as a 
function of substrate resistivity peaks for cells having 
resistivities of approximately 0.1 ohm-cm. 

The major technical findings and other aspects of progress 
made during this year of grant support are presented in the 
following sequence. 
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In Section II, we summarize our theoretical anticipa- 
tions about which of the multitude of possible physical 
mechanisms occurring in silicon junction devices predominate 
in producing the discrepancy here between theory and 
experiment. The detailed study from which these conclusions 
are drawn- appears in Appendix A. 

This study suggests that two mechanisms are likely to 
be of main importance: (1) recombination losses asso.ciated 

with defects introduced during junction formation; and (2) 
inhomogeneity of defects and impurities across the area of 
the cell. These mechanisms were among those suggested 
initially by Sah as likely candidates for study under this 
grant. In Section III, we describe the various diode test 
structures we have designed and fabricated and the measurement 
configurations we have constructed to characterize defect- 
properties and areal inhomogeneity. In addition, we report 
preliminary experimental findings. 

Section IV describes future plans for the research. 
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II . RECOMBINATION LOSSES AND AREAL INHOMOGENEITY 


As was noted in the previous section, the expectation 
of good performance has prompted interest in solar cells 
with highly-doped substrates. As an example of these 
expectations, consider from a physical standpoint the open- 
circuit voltage. Open-circuited conditions imply no current 
at the terminals. To satisfy this constraint, an opposing 
component of current must arise to balance the photogenerated 
component. In the traditional theory, this balancing (dark) 
current derives from a lowering of the potential barrier at 
the junction, which allows excess charge to flow into the 
substrate (base) where it accumulates sufficiently to sustain 
the dark current. To the first order, the potential barrier 
lowers by an amount equal to the open-circuit voltage. The 
more highly doped is the substrate, the greater is the 
barrier height at equilibrixmi; and so the more it must lower 
to produce the dark current. Thus arises the expectation 
that open-circuit voltage will increase with increased doping 
of the substrate. 

Contrary to this expectation, however, is the increased 
significance of regions near the surface that accompanies 
greater doping in the substrate. This arises for the following 
reason. The greater the substrate doping, the greater is the 
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lowering of the potential harrier. Consequently, the 
greater is the proportion of the excess charge accumulating 
in the junction space- charge region and in the quasi-neutral 
diffused layer (emitter) relative to that accumulating in 
the base. 

Thus, as substrate doping increases, the physical 
mechanisms occurring in the emitter and in the space-charge 
region gain significance. Many different mechanisms, 
traditionally ignored, occur in these regions of high doping 
and high impurity gradients. That the traditional theory 
fails to agree with experiment suggests these mechanisms 
need study. 

A major problem in understanding the operation of 
highly-doped cells lies in determining which of the mechanisms 
are dominant and which may be neglected. A second problem, 
relating to design, lies in controlling both the dominance 
and the magnitude of the phenomena via controlling the device 
structure and the steps used in fabrication, particularly 
during junction formation. 

The issue of dominance among physical mechanisms is the 
subject of the theoretical study reported in Appendix A. 

Here, in the present section, we summarize the main, conclusions 
deriving from this study. We dispense with most references 
in this section; the references for Appendix A are listed 
separately from those for the rest of the report, allowing 
Appendix A to be read as a self-contained unit. App^endix A 
deals wholly with a one- dimensional model of the solar cell; 
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in the present section, we suggest the significance of 

% 

another mechanism whose consideration produces a three- • 
dimensional model of cell behavior. 

UNDERSTANDING REACHED ABOUT DOMINANT MECHANISMS : To examine 

the issue of dominance among the high-doping mechanisms, we 
have divided them into two broad categories : 

1. .Gap shrinkage , as produced, for example, 
by band tailing, impurity-band widening 
and impurity misfit; and 

2. Altered Interband transition (recombination) 
rates arising from Auger- impact or Shockley- 
Read-Hall (SRH) processes or from electronic 
tunneling via defects . 

Which of these mechanisms predominates depends, in general, 
on the physical make-up of the device, on environmental 
conditions such as temperature, and on the aspect of cell 
performance of interest. 

In this generality we have given only qualitative 
indications about the dominance of the mechanisms . To 
provide a quantitative illustration, we have taken a concrete 
example: a phosphorous-dif fused n+p cell, junction depth 

0 . 25 . microns , impurity grade constant 10 atoms/cm^, sub- 
strate resistivity 0.1 obrnrcm. Further our attention has 
centered on the measured open- circuit voltage at 300°K. 

To analyze this device, we have extended the traditional 
analytical theory of silicon solar cells to enable inclusion 
of the high doping mechanisms. Of these mechanisms, we have 
concluded that gap shrinkage, taken alone in a one-dimensional 
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model, falls far short of explaining the measured open- 
circuit voltage. To fit the data, a gap shrinkage of 0.23 eV 
would be required for impurity concentrations only slightly 
higher- than 10 cm , which compares to our upper-bound 
estimate of 0.07 eV for such concentrations. From a physical 
standpoint, we predict gap shrinkage to be small because 
minority carriers can exist in sizable numbers in the dark 
cell only where the doping is relatively small. This is one 
of the main consequences of the piecewise linear- junction 
approximation' used in our analysis, as contrasted with the 
step- junction model customarily used in the analysis of solar 
cells . 

Of all the one-dimensional mechanisms considered,, we 
have proposed the recombination processes associated with 
defects created during junction formation to be the most 
likely candidates to explain the open- circuit voltage 
experimentally observed. From a theoretical argument based 
on the law of mass action combined with neutrality, we find 
a first approximation for the spatial dependence of the 
defects near the surface: that the defect density rises in 

proportion as the rth power of the shallow- impurity density. 
To bring the open-circuit voltage predicted by our theory 
into agreement with that observed experimentally requires 
that the exponent r take values between below 2 to about 4, 
the uncertainty resulting from certain uncertainties about : 
structure and processes that are discussed in Appendix A. 

The value r=2 corresponds to di-vacancies or other simple 
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configurations of the defects; higher values of r imply 
more complex configurations. 

Our ultimate focus in Appendix A on the spatial 
dependence of the defect density came through elimination 
of the dominance of the other mechanisms named previously, 
including gap shrinkage; nonetheless, we anticipate another 
mechanism may rule, at least in part, many aspects of 
behavior of low- resistivity cells. 

We have previously considered a one- dimensional model 
of the cell, the only coordinate of interest having been 
that measuring the distance from the surface. But the 
solar cell is a large-area device, and inhomogeneities across 
this area could play a significant role in governing the 
performance. In particular, we note the existence of a- 
spatial distribution of impurity clusters , thermodynamically 
Stable, occurring in the diffused layer [8] '. This inhomogeneity 
may result, in part, from inescapable statistical fluctuations 
and, in part, from man-made contributions that could be 
minimized. 

Combined with the one- dimensional model advance'd in 
Appendix A, consideration of the effects of this areal inhomo- 
geneity gives a three-dimensional characterization of the 
behavior of solar cells. Thus one can systematically view 
the influence of the one-dimensional mechanisms discussed in 
Appendix A in the context of inhomogeneity across the area 
of the cell. We stress that the overall effect on behavior 
is not a simple average over the area of the cell. Rather, 
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certain of the mechanisms are much emphasized, for example, 
in a region of large- impurity concentration due to 
clustering. An example is gap shrinkage, which enters as 
an exponent in detemining cell properties in a localized 
region; another example is the associated increased defect 
density and the increased recombination rates corresponding 
to it. Viewing the overall solar cell as a collection of 
sub-cells roughly in parallel one with another, we can see 
that those sub-cells with relatively high doping and defect 
density can severely degrade the performance of the device. 
Hence the areal~ inhomogeneity mechanism could play a dominant 
role and establish a basic limitation on the performance 

r 

obtainable.' 

For the cell under study, therefore, we anticipate on 
theoretical grounds that two mechanisms may likely combine 
to yield the observed open-circuit voltage:, (p increased 
recombination rates near the surface; and (2) areal inhomo- 
geneity. The experimental studies discussed in the following 
section are designed to assess these theoretical anticipations. 
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Ill . EXPERIMENTAL STUDIES 


The intent in the experimental phase of the program 
is to test, and guide the theoretical anticipations . In 
particular, the immediate concern is whether the dis- 
crepancy .between measured and predicted open- circuit 
voltage in 0,1 ohm-cm cells arises mainly from the two 
contributors singled out in Section II: areal inhomogeneity [8] 

and recombination losses, especially those introduced during 
junction formation [9] . 

The general plan of the experimentation involves 
several parts. To study the recombination losses, we 
fabricate n-on-p and p-on-n diodes , made with small area 
to de- emphasize the contribution of areal inhomogeneity on 
device performance. The current-voltage characteristics 
are measured with the effects of surface channels minimized 
by the application of voltage to an MOS guard ring.- These 
characteristics are studied for their relation, on the one' 
hand, to the open-circuit voltage, and, on the other hand, 
to the parameters describing the defect centers. The open- 
circuit voltage is measured by use of a special mask designed 
to let light penetrate a major portion of the surface area 
of the diode. The parameters of the defects (including the 
densities, activation energies, thermal emission rates and 
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their spatial dependencies) are measured by capacitance-- 
transient and related methods due to Sah and his co- 
workers [9-12]. By applying these methods to Schottky-; 
barrier diodes of appropriate structure, one may characterize 
the defects present before junction formation, and, hence, 
infer the defects and the defect properties deriving from, 
the processing used in pn junction formation. 

To study areal inhomogeneity, we routinely map the 

current- voltage characteristics over the area of wafers that 

are processed. This is done by making measurements on small- 
-2 2 

area (~ 10 . cm ) diodes scribed from a large-area (1.25” 
diameter) wafer . On these small-area devices, one can also 
measure .the defect properties and the photo-voltage, and 
thus secure an areal mapping of many parameters that in- 
fluence -cell behavior. Alternatively, one can map the photo- 
voltage by probing before the wafer is scribed [13] , 

Two aspects of the processing are noteworthy. First, 
for surface preparation, we have built a tilted, teflon 
beaker, which rotates by a speed-controlled motor. The 
beaker enables a chemical polishing of silicon wafers that 
reduces the damage present on the surface. Second, we have 
devised processing,, including diffusions, for fabri-' 
cation with temperatures never exceeding 900 “C. The intent 
here is to study devices with low defect densities,, made by . 
processing similar to that employed in conventional solar- 
cell fabrication. 
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Preliminary experimental findings on both np and pn- 
devices showed good characteristics of the small-area 
diodes made on 0,1 ohm-cm substrates. Saturation current 
densities of about 10 Amper.es/cm and slope factors of 
about 1.1 .were commonly seen. The experimental configurations 
for the transient-capacitance measurements were burlt and 
shown to work-. By use of a tungsten lamp, we secured an 
indication of the open- circuit voltage on the small-area . 
structures; 590' millivolts was not uncommon, which. correlated 
well with predictions based on the observed dark-current 
characteristics joined with an estimate of the short-circuit 
current. The mappings revealed considerable areal inhomogeneity 
of the parameters describing the dark current and of the open- 
circuit voltage. This suggests the presence of a man-made 
contribution to the inhomogeneity. 

The effort spent during this first year of-grant support 
has put us in a good position for systematic and detailed 
experimental studies during the year of renewed support. 
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IV . FUTURE PLANS 


During this period of support, we spent much effort' 
trying to sort, on theoretical grounds, the prime contribu- 
tors to the low experimental open-circuit voltage from' 
among the multitude of different phenomena that occur in 
silicon devices. Section II sets forth the conclusions 
deriving from this effort: areal inhomogeneity and re- 

combination losses near the surface are the likely prime 
contributors for the cell under study. 

In the immediate future, our efforts will concentrate 
on an experimental assessment of these theoretical anticipa- 
^tions. As is implied in Appendix A, to enable such an ■ 
assessment requires first filling the present gaps in under- 
standing about the recombination losses in substrates of 0.1 
ohm-cm resistivity, including the total absence of informa- 
tion about the spatial dependence of the lifetime. Hence' 
our initial .attention will focus on characterizing the defect 
properties in the substrate produced by the processing used 
in pn junction formation. For this study, the decreasing 
sensitivity of the transient-capacitance methods with increased 
substrate doping may pose experimental problems, which will 
require special experimental strategies . ■ • 

In the farther future, still concentrating on cells with 
0.1 ohm-cm substrate resistivity, we plan to examine further 
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the role of the recombination losses and other phenomena-- 
occurring in the highly- doped emitter region. Throughout, 
we intend td monitor the areal inhomogeneity. Open- circuit 
voltage will continue to be the main concern. 

In the far future, we foresee an engineering phase of 
this study. This involves a systematic linking of ’the 
defect properties, the areal inhomogeneity, and whatever 
other physical mechanisms may emerge significant to the 
processing used in fabrication. The intent is that the 
scientific understanding gained shall underlie the develop- 
ment of batter solar-cell performance. 
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APPENDIX- A: DETAILED CONSIDERATION OF THE ONE-'DIMENSIONAL 

PHYSICAL MECHANISMS* 

INTRODUCTION 

Recent work [l,2] based on the traditional analysis of the 
silicon junction solar cell predicts an improvement in perfor- 
mance as the substrate doping is increased. But, as noted by 
Brandhorst [l], experiment failes to agree with theory. For 
example, measured open-circuit voltage peaks at a substrate resis- 
tivity of approximately 0.1 ohm-cm, and shows a maximum value of 
610 mV for this resistivity [3], about ten to fifteen per cent 
below the predicted value [3j. 

Two possible sources of such disagreements are apparent. 

First, they may originate in the approximations made in tradi- 
tional analysis; and niimerical solution by computer can help 
identify this kind of difficulty and remedy it. This is the 
approach and main focus of the recent work of Fossum [4] and of 
Dunbar and Hauser [5,6]. Second, hovrever, the disagreements may 
partly owe their origin instead to the neglect in traditional 
analysis of certain physical mechanisms accompanying high doping. 
This possibility is the subject here. 

Our intent here is to describe mechanisms ignored in the 
customary view of cell operation, and to indicate the degree to 
which they can influence solar-cell characteristics. These 


This Appendix derives from part of the contents of the paper: 

F. A. Lindholm, et al . , Record of Eleventh IEEE Photovoltaic 
Specialists Conference, pp. 3-12, (1975). Though the general 
conclusions of the two works agree, the Appendix contains several 
noteworthy additions and alterations. 
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mechanisms all arise from high doping; and some can substantively 
determine the kinetics of holes and electrons in the highly-doped 
diffused layer and in the junction space-charge region. For the 

t 

most part, the physical mechanisms to be discussed tend to degrade 
the performance of solar cells. Hence they impose fundamental 
limitations on the performance obtainable. 

As background, for our treatment, we begin with a generaliza- 
tion of the customary view of solar-cell operation , with emphasis 
given to the origin of the open-circuit voltage. Next we propose 
the additional mechanisms and comment qualitatively on their rela- 
tion to temperature and to the physical make-up of the cell. Then 
follows an extension of the traditional analysis and a comparison 
to experiment. Thus, for a one-dimensional model, we indicate 
which of the additional mechanisms dominate in determining the 
open-circuit voltage at 300°K of a cell with 0.1 ohm-cm substrate 
resistivity. 

A GENERALIZED VIEW OF TRADITIONAL ANALYSIS 

In the traditional method of analysis, as illustrated in 
Figure A-1, the cell is sectioned into three regions: the quasi- 

neutral portion of the substrate (or the base) , the quasi— neutral 
portion of the diffused layer (or the emitter) , and the junction 
space-charge region (or the SCR) . Further, one regards the cur- 
rent density at the terminals as the superposition of two oppositely 
directed components, 

j - jj, - (1) 
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Fig. A-1 


Qualitative sketch of mobile charges in dark cell, suggesting 
dominance of Qg, as :n traditional theory. Charge Qg anticipates 
one high-doping effect. Charge is not shown. ’ 



in which denotes the total spectral current density generated 
by the solar illumination, and denotes the dark current density 
that V 70 uld occur in the absence of illumination. Traditional 
analysis, as reviewed by Brandhorst [l] and Stirn [7], then bases 
the determination of on a step- junction model and on simple 

diffusion theory. Here we avoid these idealizations by a more 
general view which, as is shown later, enables the inclusion of 
high- doping me chan i sms . 


We regard as the sum of three components: 

+ '^SCR ^ ''e • 


(2) 


VJe think of each component as deriving from the recombinat ion-diffusion- 
drift kinetics of the mobile carriers in the three regions of the device,' 
and thus regard each as the ratio of an excess mobile-carrier charge to a 
time constant pertaining to the region under consideration. .Hence,, 


•^B = 

QbAb ' 

(3.) 

■^SCR 

” ^scr'^^scr ' 

(4) 

= 


(5) 


Figure A-1 conveys this view of the origin of the dark current. 
Depending on the process dominating the kinetics of the mobile 
carriers, the time constant will often have the physical interpreta- 
tion either of a lifetime or of a drift-diffusion (transit) time or 
of. a combination of the two. The view expressed in Eqns . (3)- (5) is 

general, for a one~dimensional model, inde'pendent of the doping 
profiles of the mechanisms controlling cell behavior, , 

Traditional analysis leads to the conclusion that Jg constitutes 
nearly all of the dark current density at 300°K; that is, for 
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open-circuited conditions 


<«> 

Figure A-1 pictures the reason for this. The thickness of the 
base region by far exceeds that of the emitter or the junction 
space-charge regions, and consequently the charge Qg apparently 
much exceeds. Q„ or Hence the base current will tend to 

ht bUK 

predominate provided one supposes, as is customary, practically 
equal lifetimes in the three regions and a sufficiently large 
drift time for the emitter. 

The traditional analysis* then models the dark current density 
using Shockley's diode theory [8]: 

Jg = jQ[exp{qV/kT) - l] = Jq exp (qV/kT) (7)’ 

Here denotes the saturation current density, which is propor- 
tional to the diffusion length Lg and is inversely proportional - 
to the doping, of the base : 






From Eq. (7) and (8)- it follows that 


Qg exp (qV/kT) 


( 8 ) 


(9) 


Shockley's theory [8] applies to an idealized p.-n junction.' Actual 
junctions show a current deriving from net recombination in the 
quasi-neutral regions that varies in proportion as exp(qV/mkT), 
in which m often is as much as ten or twenty per cent greater 
than unity. Inclusion of m?^l in the arguments presented here 
causes no substantive change, however; and thus, for simplicity 
of discussion, we assume m=l throughout. 
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which is the charge of the minority carriers in the base appeiaring 
irt Eq. (3) and illustrated in Figure 1. Iii the notation used 
above, V stands for the applied potential across the junction space- 
charge region, and the other symbols have their usual meanings. 

This conception of solar-cell operation enables an easy inter- 
pretation of the open-circuit voltage . To conform to open- 
circuited conditions, the minority charge in the base Q_ must 
rise until is large enough to exactly balance This requires 

the applied component of the barrier potential to rise until it 
reaches 

Vqc s kT/q (10) 

in the forward direction. Examination of Figure A-1 and of 
Eqs, (8) - (10) shows that two factors compete to control 
as the substrate doping Ng increases. On the one hand, an increasd 
in tends to depress the boundary value of the minority carriers, 
and thus Vqq must increase to let more minority carriers flow into . 
thfe base. On the other hand, empirical evidence exists that the 
lifetime drops [9] as doping increases, which demands less Q_ to 
constitute a given J_; this tends to decrease . In the tradi- 
tional theory, however, the first tendency dominates and is 

predicted to increase as doping rises [],2]. 


22 



Having established this background, we can now pinpoint flaws 
in the traditional method of analysis that may help explain its 
large predicted values of for highly-doped cells. First,. . 

recall that and were neglected in comparison with J^. Their 
presence in significant values could by itself explain the disagreement 
between experiment and theory. This one can see -mathematically 
from Eq. (10).; Or, one can view this from a more physical stand- 
point, using Figure A-1. Note that the presence of significant 
excess currents or and hence of significant values of 


"E 


and QgQ, diminishes the value of needed to balance J,. . Since 


B 


according to Eq. (9) , depends exponentially on the presence 


OC' 


of excess currents will therefore depress Vq^^. 

But Figure- A-1, which is scaled for. a substrate of ordinary ' 

.doping, suggests that the excess currents are wholly negligible. 

What changes occur with high doping that give rise to sizable 

values of such change arises from the observed' 

* ^ 

drop in lifetime in the substrate [9]' as substrate doping increases. ' 
This, for example, tends to give greater importance. 

I 

But other changes occur too. As one. can see from -Eqs. (3) - (5),: 
any mechanisms accompanying high doping that increase or' Q„ 

(as suggested by the dashed shape in Figure A-1,) or that decrease 
the corresponding time constants can produce sizable excess cur-- 
rents, and can contribute to the cell performance seen experi- 
mentally. Many such mechanisms exist, and these are discussed 
in the following section. 


The processing associated with junction formation might also . 
affect the lifetime in the base. The experimental studies described 
in Section III address this question, among others. 
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ADDITIONAL PHYSICAL' MECHANISMS 

The traditional analysis can be viewed from another vantage 
point, which enables a concise s\iminary of the physical mechanisms 
it contains. As we have noted, the traditional analysis uses 
Shockley's theory of the pn diode [8]. Hence, fundamentally, it 


relies on the equations underlying Shockley's analysis. 

which are: 



(11) 

or 


(11a) 


Jp = -qDpVP + qUpEP 

(12) 

or 

Jp = -qypV(|)p ^ 

(12a) 


j = + Jp 

(13) 


° ^ ^EXT “ ^ 

(14)' 


0 = - - Rj, 

(15) 


-V^4, = VE = I (H-P+N^-Npp) - 1 

(16) 

These. 

equations, joined with the auxiliary relations. 



N = n^ exp[q({f)-(j)^)/kT] 

(17) 


P = n^ exp[q(^p-(j))/kT] , 

(18) 


constitute the starting point not only for the traditional analysis 
but also for the computer studies of' cell operation now beginning’ 
to appear [4-6"] . 

The first of these equations expresses the electron current 
density as the sum of a diffusion component, proportional to 
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the gradient of the density n of mobile electrons, and of .a drift 
component, proportional to the electric field E. The constants 
of proportionality involve the mobility and the dif fusitivity 
D^, together with the magnitude of the electron charge q. In 
Eq, (11a) appears an equivalent way of expressing the electron 
current density, in terras of the gradient of the quasi-Fermi, 
level <|)j^ of electrons. Eqs.' (12) and (12a) are the counterparts 
for the mobile holes; and Eq. (13) expresses the total current 
density as the sum, of electron and hole components. By Eqs. (14) 
and (15), the continuity equations, one acknowledges that the time 
rate of change of mobile carriers within an infinitesimal (but 
macroscopic) volume arises from two sources: from the rate of 

flow into the volume and from the net rate (G„^,„ - R) of generation 
within the volume. Eq. (16) is the Poisson equation that deter- 
mines the electric field in terms of the charge density. The 
components of the charge density include the densities, and 

of the acceptor and donor impurity ions. The dielectric permit- 
tivity e enters as usual. 

In the analysis of solar cells, the term describing 

external generation receives special attention; it gives rise to 
the total spectral current density J generated by the illumina- 
tion. The recombination rates and Rp are customarily modeled 
eithe'r simply in terms of carrier lifetime of by assuming a single 
SRH defect center lying near mid-gap [10,11]. As we have noted, 
the traditional analysis treats the lifetime as spatially inde- 
pendent, though its dependence on doping in the substrate [9] 
has recently been included [3]. 
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Eqs. .(11) - (18) constitute the physical mechanisms conven- 
tionally included in the analysis of solar-cell operation. By 

4 

computer solution, one can take account of the effects of these 
mechanisms on cell operation with near exactness. The method of 
traditional analysis obscures their effects through the ideal- - 
izations made, though it can give useful physical insight. 

We now discuss additional physical mechanisms, not customarily 
included, which can occur in highly-doped cells and can signifi- 
cantly influence cell characteristics. As we shall see, these 
mechanisms generally tend to increase the excess currents and 
“^SCR Eqs. (4) and (5) . They do so either by increasing 

the excess mobile charges, Q„ and Q„™, or by decreasing the 

hi oCK 

related time constants, and' 

The mechanisms to be discussed fall into two broad categories: 

(a) Distortion of the energy-band edges and production thereby 

of effective gap shrinkage (or an effective intrinsic density) 
and quasi-electric fields . Gap shrinkage tends to increase 
the excess mobile charge, particularly Q^. The quasi-electric 
fields can affect both excess mobile charge and the transit 
times, particularly in the emitter. These mechanisms alter 
Eqs. (11), (12), (17), and (18) of the basic set underlying 

the traditional analysis of solar cells. 

(b) Alteration of the recombination rates (rates of interband 


transitions), denoted by R and R in Eqs. (14) and (15). 
n p . 

This can directly influence the time constants x„ and x,,™, 

iii oUK. 

and, through this influence, the excess mobile charges Q^, 
®SCR* 
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We now discuss these two categories of mechanisms. 

Effective Gap Shrinkage 

At high ■ doping concentrations , various deviations from the 
conventional band structure can occur, such as: 

(a) The band edges need no longer be sharp; rather, states can 
extend beyond the band extrema describing lowly-doped sili- 
con, forming thereby so-called band tails [I2j; 

(b) The impurity states can broaden into an impurity band [13] 
with increasing impurity concentration; at high levels of 
doping these states can become delocalized; 

(c) A spatial variation of the energy gap can. result from macro- 

scopic lattice strain introduced by a high-concentration 
diffused impurity profile, the strain arising from the. 
misfit of the impurities in the silicon lattice ' [14j . ' 

As Figure A-2 illustrates qualitatively, each of these devia- 
tions tends to produce an effective shrinkage AE^ of the energy 
band gap. Because the intrinsic density n^ depends exponentially 
on the band gap, the gap' shrinkage can be imbedded in an effective 
intrinsic density [15], which then appears directly in Eqs. (17) 
and (18) of the basic set. Note that spatial dependence of the 
doping concentration implies spatial dependence of the effective 
band gap. 

To examine qualitatively the effect of gap shrinkage on 
the conduction processes described in Eqs. (11) and (12) of the 
basic set, we must distinguish between whether the states intro- 
duced by high doping are localized or delocalized. Thus, bor- 
rowing from the concepts used in models describing amorphous 
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semiconductors, we identify a critical energy, the Cohen 
mobility edge [16] , which separates localized from delocalized 
states; the mobility edge plays a role in describing charge 
transport in disordered solids similar to that which the band 
edge plays in periodic solids such as lightly-doped silicon. 

Then, for delocalized states, conduction can occur by drift 
and diffusion; if the variation in doping produces spatial 
dependence of the band gap or of the electron affinity, quasi- 
electric fields [17] can influence the transport of electrons 
and holes. For localized states, conduction can occur from 
hopping between neighboring states [18] due to thermally acti- 
vated transfer of the trapped charges in the Mott model [19] 
or by tunneling in the Anderson model [20]. 

For ordinary temperatures, only the delocalized states will 
contribute appreciably to current flow. But, in such states, 
notice that three forces can act to produce the flow of holes 
and electrons: the statistical force due to diffusion and the 

Coulomb force due to the electric field — just as in lightly- 
doped silicon — plus the effect of the quasi-electric fields 
(another statistical force) . Thus Eqs. (11) and (12) each 
require one additional term; however, Eqs. (11a) and (12a), in 
terms of the quasi-Fermi potentials, remain valid. 

Interband Transition Rates (Recombination Rates) 

Transitions between localized and delocalized states that 
constitute recombination-generation-trapping processes occur 
via several energy-momentum conservation mechanisms , including 
phonon-assisted , photon-assisted , and Auger-impact processes [21]. 
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The rates of such transitions depend on the density of the 
defect centers, and the corresponding energy levels, capture 
.cross sections and emission rate constants. By a combination of 
experimental methods [22], one can establish values for all of 
these in the junction space-charge region and the base of the 
cell. In the highly-doped emitter, one expects the defect density 
to rise sharply near the degenerately-doped surface. This expec- 
tation arises from the strong' suggestion of experiment [23]^ that 
defects in silicon device processing are generated at the strained 
disordered surface layer, from where they diffuse rapidly into 
the bulk. These experiments further suggest that the defects 
are silicon vacancies or vacancy complexes. Thus, using the 
analogy to chemical reactions reviewed by Fuller [24], one finds 
the following first-order model for the dependence of the &efect 
density ^TT on the doping in the- emitter and space-charge 
regions . 

N^^(x) = k[n^j^(x) + (19) 

which applies for sufficiently large 

In deriving this model, one assumes approximate thermal 
equilibrium during the diffusion of impurities; this justifies 
the law of mass action, which one combines with the neutrality 
condition. The different values of the exponent r derive from 
different hypotheses about the dominant chemical reaction. For 
reactions corresponding to the formation of di-vacancies, for 
example, the choice r = 2 provides an approximate fit. 

Aside from the processes mentioned one other kind exists 
that can influence the effective interband transition rates. 
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This kind involves tunneling via deep-level defect centers . ■ 

Among the many processes of this general type [25]/ we shall 

describe here a particular two-step electronic process. This 

process assumes an electron i-n the conduction band is first • 

captured by a localized deep level and subsequently tunnels 

to the valence band or to the. valence-band tail states;. Sah has 

¥ 

demonstrated this as the excess-current mechanism in gold-doped 
silicon' tunnel diodes [25]. 

QUALITATIVE REMARKS ABOUT WHICH MECHANISM DOMINATES 

Consider first the physical mechanisms that affect the 
interband transition rates- In the quasi-neutral substrate (the 
base) of the cell, description of these rates by a simple life- ‘ 
time, assumed spatially independent, may suffice as a first- 
order approximation, although experimental evidence is wholly 
lacking at present. We make this assumption for simplicity. 

In the junction space-charge region, the simplest model, due to 
Sah, Noyce and. Shockley [26j, assumes that electronic transitions 
between the conduction and valence band occur via localized 
states of a single energy level lying deep within the forbidden 
band. These two processes produce the base current and the SCR 
current of the traditional analysis, as defined in Eqs. (3) and 
(4) . 

■ By. .physical, reasoning, one can see what conditions tend to 
yield dominance of one of these currents oyer the other. The 
base current results from carriers that have surmounted the 
entire barrier of the junction region. To contribute to the SCR cur- 
rent, however, the carriers need to have climbed only part of 
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the barrier. Therefore any conditions that increase the barrier 

* h 

height or decrease the climbing energy of the carriers tend to 
cause dominance of SCR current. Hence the dominance of SCR cur- 
rent is more likely for loads approximating short-circuit rather 
than open-circuit conditions and for lower rather than higher 
temperatures. For very' low temperatures, one may see the appear- 
ance of the processes [25] that involve tunneling, because tun- 
neling enables holes and electrons to reach nearly the same 
spatial position v;ithout needing to climb the barrier. A steep 
junction also favors tunneling processes. 

In the quasi-neutral diffused layer (emitter) , Auger proces- 
ses may compete significantly with Shockley-Read-Hall (SRH) proces 
ses near the surface, where majority carriers exist in great • 
niombers. Near the edge of the junction space-charge region,- 
where doping is relatively low, our calculations indicate that - 
SRH processes will dominate band-to-band Auger processes-. For 
the present, the relative role of Auger recombination via defect 
centers cannot be assessed because of the absence of experimental 
studies that would fix the values of the. pertinent rate constants. 

Consider now effective gap shrinkage and its characterize- , 
tion. Recently strides have been made toward including the 
effects of gap shrinkage (often characterized by an effective 
intrinsic density) in explaining the operation of a bipolar 
transistor [27]. • From this inclusion has come a closer agree- 
ment between experiment and theoretical prediction than was 
obtained previously, particularly for the forward-active current 
gain and the gain-bandwidth product. This suggests direct 
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applicability to solar cells of the model for gap shrinkage .used 
in transistor theory; but uncertainties exist that call for 
Caution. 

The model for gap shrinkage now used in bipolar transistor 
theory is based/ essentially, on the Kleppinger-Lindholm extrap- 
olation [28] of the characterization- of Morgan [l3] for impurity- 
band widening and of Bonch-Bruyevich [12] for band tailing. No 
direct .experimental evidence has yet appeared to assess the degree 
of fit for silicon of the theories of Morgan and Bonch-Bruyevich; 

^ ft 

and Kleppinger and Lindholm intended their extrapolation to serve 
only as a first-order approximation. Further, the usage put to 
this extrapolation in the computation of transistor behavior has 
evolved quite different dependencies of the effective intrinsic 
density (or effective gap shrinkage)' on the doping density. This 
one can see by comparing dependencies appearing in. the various 
treatments [27], Finally, the theory treats the states in the 
Kleppinger-Lindholm extrapolation all as delocalized states.' 

This adds to the exaggeration of gap shrinkage inherent already 
in the extrapolation. 

Thus first-order models for gap shrinkage exist that prob- 
ably could bring the predicted and measured performance of cells 
into better accord. But we conclude that the uncertainties pre- 
sent in the existing theories . rule 'out for the present a defini-.- 
tive inclusion of the effects of' this mechanism. In assessing 
the relative significance of gap shrinkage, therefore, we shall 
try to evade the uncertainties of- the existing theories by 
relying on simple estimates.. This ppssibility we explore in the 
next section," 
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EXTENSIONS OF THE TRADITIONiu:. ANALYSIS 

Here we extend the method of the traditional analysis to 

estimate the effects of spatially-dependent gap shrinkage and 

recombination rates. We do this for a cell of 0.1 ohifi-cm res is-' 

tivity operating at a temperature of 300°K. For such a cell' at 

this temperature, consistent with our discussion in the fore- • 

going section, we anticipate that gap shrinkage and altered 

» 

recombination will be the significant additional mechanisms 
acting to determine cell performance. Our interest in a 0.1 
ohm-cm cell arises from the maximum it shows in measured open- 
circuit voltage [3]. For concreteness, we assume a phosphorous- 
diffused n+p cell with a surface concentration of 10^®cm“^ and a 
junction depth of 0.25 microns. 

As in the traditional analysis, we regard the current at 
the terminals as the sum of and in accord with Eq. '(1) . 
Because gap shrinkage occurs appreciably in the emitt^ region 
only, which constitutes a small fraction of the volume of the 
entire cell, we assume that gap shrinkage negligibly affects J- . 

li 

Thus we use the traditional value for J- [9] applicable to semi- 

1j 

infinite bulk material. 

As in Eq. (2), we regard as the sum of three components, 
and we now seek to calculate these. According to Eqs. (3) - (5), 
this involves calculation of the' excess mobile charge and the 
time constant characterizing each of the three regions of the 
device. 
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The Base Region 

Treatment of the substrate (base) region requires no changes. 
' from traditional analysis. Consistent with experimental data on 
• diffused devices [23], we assume a -spatially-independent defect 
density and lifetime in the region; and, consistent with 
■ theories of gap shrinkage [28], we assume negligible shrinkage 
for a doping concentration, 5(10^^) cm“^, v/hich corresponds 

to p-type resistivity of 0.1 ohm-cm. Thus, 

°E ^ [qJ'Bfaj/MftA e^PlqVAT)] 

E ■ ■'b 

Substitution of - 10 microns, which is " consistent with recent 

o 

data [ 9 ] for diffusion length Lg against doping,* and combination 
with Bq. (10) then yields 

Vqc - 680 - mV - . .(21) 

This agrees satisfactorily with the predictions of traditional 
; theory [3]*; the maximum measured value is 610 mV [3.]. 


There is scatter in these data of about one order of magnitude 
for a resistivity of 0.1 ohm-cm. Use of 100 microns in- this 
calculation- would yield an open- circuit voltage of about 740 .mV', 
which would not change qualitatively the conclusions to» be drawn. 
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Junction Model 


Now consider the junction space-charge region. At open- 
circuit conditions the junction barrier is forward biased by a 
voltage approximating According to depletion theory, the 

width of the junction space-charge region narrows to a value 

lower than its width at thermal equilibrium. As a first-order 
approximation in calculating we shall assume a linear ■ 

junction space-charge region; that ^ is. 


^Dd“ ^AA= -a(x-Xj) (22) 

For junctions formed by the conventional processing used in 

transistor fabrication, the doping gradient "a" does not exceed 
23-4 

about 10 cm ,’ according to recent measurements on shallow 
phosphorous-doped emitter junctions [29] . For junctions formed 
by the lower- temperature processing used in solar-cell fabrica- 
tion, we assume this value constitutes a reasonable upper bound 
for the doping gradient. 

In bipolar transistor theory, a linear junction model [8,30] is 
often used as an approximation for the emitter junction [29,31,32], 
even for devices with very shallow and heavily doped emitters. 

This approximation serves as a better model for solar-cell 
junctions than does the step- junction approximation commonly used, 
because a step junction can severely exaggerate the influence of 
gap shrinkage on cell behavior, as the discussion to follow will . 
imply. 

The linear- junction approximation itself, however, proves 
incapable of -providing good estimates of various aspects of 
behavior; for example, it underestimates the strength of the 
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built-in field and the doping density at the emitter edge of 
'the space-charge region. Thus, where needed in the course of 
•the treatment below, we' shall indicate certain refinements of 
■this approximation. 

Extensions for the Junction Space-Charge Region 

Assuming a linear junction, we now use conventional 

calculations of W which derive from solution of Poisson '-s 

scr 

equation under the depletion approximation . For a doping gradient 
of 10^^ cm these yield at T = 300°K 
■ lOOOAO 


which decreases to about half of this value for 
V - 610 mV. 

The thinning is less than this because the presence of mobile 
carriers invalidates the depletion approximation [30] . Moreover, the 
asymmetrical departures of the actual profile from that of Eqn. 

(22) imply that less than one-half of the total thickness of the • 
space-charge region occurs on the emitter side of the. metallurgical 
junction. But we content ourselves with estimates here, recognizing 
that several Debye lengths (which mark the width of the blurs 
separating space-charge and quasi-neutral regions) occupy several 
hundred Angstroms of material for the doping levels' of interest. 

Using the estimate = 1000 A° , together with the assumption 

oLK 

23 —4 

of' a linear junction and of "a" = 10 cm , we find the maximum 
value of the doping in the space-charge region: 


Nj^(Xg) - 10^^ cm ? 


(23) 
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in which denotes the boundary separating the quasi-neutral 

emitter from the space-charge region. This constitutes a lower- 

bound estimate because the profile on the' emitter side of the 

rises more rapidly than does the linear profile of Eqn. (22) ; 

and, as shown by calculations based on a piece-wise linear 

approximation of the actual profile, the maximum density may be 

as much as twice that given in Eqn. (23) . Neither experimental 

evidence [3 3] nor theory [28] , however, suggests appreciable gap 

18 — 3 

shrinkage for doping concentrations less than 2(10 )cm ; and 
thus we neglect this mechanism in the junction space-charge 
region. This neglect, which the treatment in the next section 
substantiates in more detail, is one central consequence and 
simplification yielded by the linear- junction model and its piece- 
wise linear refinement. 

According to Eq. (19) , however, the defect density and hence 
the lifetime can vary significantly with position in the space- 
charge region. In using Eq. (19) to demonstrate this, we assume 
it to be valid also for small Nj^ ‘ This additional approximation 
simplifies the estimates to be made and allows a one-parameter fit 
to' data. From a physical standpoint, it has the drawback, however, 
of suggesting the same kind of defect predominates in the bulk as 
near the surface, which is questionable [23].' To simplify the 
estimates further, we assume the lifetime varies inversely with the 
defect density; this approximation probably introduces only small 
error in the absence of significant Auger recombination via traps. 

Thus, in the space-charge region. 
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( 24 ) 


N^^Cx) _ Tb _ ^Ndd(x), + 5(10^'^)]^ 
^TTB ^SC .5(10^^) 


This dependence is shown qualitatively in Figure A- 3 with 
(x) given by Eq. (22). The maximum value of this ratio occurs 
at the emitter edge of the spaceTcharge region, and is 3 , which 
implies about an order of magnitude decrease in lifetime for r=2 
and two orders of magnitude for r=4. A simple linear average 
across the width of the space-charge region produces the 
following values for the ratio : 2 for r=l, 3 for r=2, 

7 for.r=3; and 16 for r=4. 


To specify the current W©' identify the pertinent charge 

Qg^ from the Sah-Noyce-Shockley theory [26] and choose an appropriate 
average lifetime Tg^ for the region-: 


^SC [qn^Wg^Iexp [qV/2kT] 


For the current density Jg^^ in this cell, the main departure 
from traditional theory caused by heavy doping in a one-dimensional 
model results from the rise in defect density and fall of lifetime 
near the highly-doped surface. 

Extensions for the -Emitter Region " 

The quasi-neutral emitter region (diffused layer) of the cell 
under study contains very high doping concentrations. Nonetheless, 
we shall indicate that the resultant gap shrinkage in a one~ 
dimensional model fails to give enough deviation from traditional 
predictions to fit empirical data, whereas spatially-dependent 
lifetime does suffice. This suggests that, for the cell under 
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cpnsideration, at a temperature of 300°K, the main additional 
mechanism in a one~dimensional model governing J„ of the emitter 

hi* ' 

is spatially-dependent lifetime. 

Our approach here will involve estimating upper bounds bn 
the gap shrinkage possible. The upper bounds to be used will 
depend strongly on the spatial distribution of the holes in the 
emitter, and our treatment will start with that issue. 

The key insight lies in recognizing that the built-in 
electric field produced by the doping gradient in the emitter 
packs the minority carriers (holes) very near the .space-charge- 
region edge. Kennedy and Murley [31] were the first to identify 
this compression as a main contributor to the injection efficiency 
of bipolar transistors. 

Their calculations are based on a refinement of the linear 
junction model whose power-series approximation to the profile of 
the emitter side of the metallurgical junction is meant to provide 
a good estimate of the built-in electric field in the quasi-neutral 
emitter. The calculations show the holes confined close to the 
space-charge edge x , within a region ten percent of the thickness 
of the entire quasi-neutral emitter. Specifically, they. show that 
within this region the hole density drops by an order of magnitude 
from the value p(x^) at the edge of the space-charge region. It 
drops with the spatial distribution shown qualitatively in Figure 
A-4, which implies that practically all holes vanish within .a 
region thirty percent of the emitter thickness. 

Although the calculations apply to a device different from 
the cell under study, the differences further support the contention 
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Fig. A-4 Approximate distribution of minority carriers in quasi-neutrai 
diffused layer (emitter) of dark cell (assuming the quasi- 
fields are negligible near where Nj^ is relatively small;) . 
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that the holes, exist in appreciable numbers only in a thin region, 
.near x^„ Kennedy and Murley assumed a junction depth = 1.75 
microns, considerably deeper than the value Xj = 0.25 microns 
assumed here.. But a shallower junction and consequent higher 
built-in field would only tend to compress the holes more tightly. 
Kennedy and Murley assumed a spatially-independent lifetime, and 
a lifetime that varies' in accord with Eq. (19) would again only 
add to the compacting. 

Thus we conclude that the holes in the emitter of the cell 
occupy at most ten percent (200 A°) of the quasi-neutral emitter, 
and are confined to the lowest-doped region of the emitter. In 
all regions of the dark cell significant to its performance, 
therefore, gap shrinkage Is small. 

To place an upper bound on the 'effects of gap shrinkage, we 
assume the rectangular distribution .shown dashed in Figure Ay4. 
This distribution conflicts with the demands of hole diffusion, 
but it is meant only to clarify which mechanism dominates rather 
than .as a tool for decisive calculation. For the thickness of 
the base of the rectangle, we choose = 200A° consistent with 
the discussion above and recognizing that this choice corresponds 
to about two Debye, lengths for this doping; hence choice of a 
smaller value obscures physical meaning. 

In this region of hole occupancy in the emitter, if we 

continue to assume the linear- junction model of Eg. (22) , we find 

18 "^3 

that N will have ri'sen to 1,2(10 ) cm at x*^ which defines 

UJJ ll/ 

the edge of the hole distribution nearest the surface. As 
discussed in the foregoing section, ,a piece-wise linear approxima- 
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tion to the doping profile implies =; 2(10 ) /cm , which 

we choose for subsequent calculations . 

Estimation of from Eq. (5) requires estimation of and 
Tg. For Qg, we .extend conventional theory by inserting in the 
customary expression for nj^(x^) a band gap altered by the 
shrinkage AEg(Xg) . Thus, 

Qg = {qLgn?(x^)/[Nj^(Xg) - N^] 


T 2 

qhEUi 


Nd (Xg) . 


N 


^AEG(x^)/kT ^qV/kT 


A 


(26) 


To give an upper bound for AEg(Xg), note the empirical evidence 

of Pearson and Bardeen [33] that the ionization energy of shallow 

18 

dopants in silicon vanishes for concentrations between 10 and 
19 -3 

10 cm . As a rough estimate, we interpret this to imply that 
the conduction band extends downwards by the ionization energy of 
phosphorous: - AE^ = 0.044 eV. We add this with Fritzche’s 
description of band-tail states [34] — giving a shift upwards of 
the valence band of AE^ - 0.025 eV. This yields AEq(x^} O'. 07 eV 
2.. 7kT and exp (AE'^ (x^ /kT) - 15 . • 

To estimate t^Cx^), we use the reasoning of the foregoing 
subsection, adding the consideration that the assumed shift AE^ 
in the effective valence band- increases the hole capture cross 
section by about a factor, exp(AE^/kT) =; e. Thus, 

e(5)^ . (27) 

which gives: 14 for r=l; 68 for '^=2; 340 for r=3; 1700 for r=4 . 

This departure from traditional theory compares to- the upper- 
bound estimate of 15 for the departure arising from gap shrinkage. 
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In accord with Eq.. (5) 


~ [NnCXg) " N^]t, 


A" '•B 


qV/kT rr / 1 

X [Lg/Xg] 


AE (xM 

X [,exp ] X [e{5) ] 




( 28 ) 


The first term is the for a shallow- junction cell that 

results from application of traditional theory to an abrupt- junction 
model. . The next three terms are correction factors arising .from 
three different physical mechanisms; The first of these F, whose 
value is approximately F = 0.1, accounts for the compacting of 
the hole density. by the built-in electric field. The second G, 
which has an upper bound G=15, accounts for gap shrinkage. The 
last S accounts for a spatially-dependent lifetime, and has the 
values listed above for different values of r. 

Implications of Measured 

From a physical point of view, the open-circuit voltage 
arises to establish just large enough to balance Of the 

three components of Jj^, will be highly significant at lower 

temperatures, for the reasons noted earlier. For T = 300®K, 
calculations based on our previously stated assumptions show it 


to be less than J,, for V - V,.,„ . 

j5 UU 


Thus we ignore for the 


present discussion. 

From Eq. (10) , therefore, to fit empirical data. 


■^OE ■*' " exp(qAV^^/kT) J, 


OC' 


OB 


(29) 


implying 


Jog 14J^g 


(30) 
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Here denotes the difference between the measured value of 

•610 mV and the theoretical value of 680 mV calculated from 
traditional theory assuming J_ >> J_. and J_„ denote the 

saturation currents for J„ and J„. 

ill 

Combining Eqs. (20) and (28) gives: 

where 

^OET^^OB = ( 32 ) 

is the ratio predicted by traditional theory and FGS represents 
the correction factor arising from the three physical mechanisms 
described earlier. To fit empirical data, FGS - 700. 

The factor F, representing the effect of the built-in 
electric field in the emitter, is approximately 
F = L^/Kj. - 0.1 

The factor G, which accounts for the effects of gap shrinkage, has 
an upper bound of 15. To fit the data, therefore, the factor S, 
which describes the influence of spatially-dependent lifetime, 
must be at least 466. This implies an exponent r=3 in Eq. (19) 
describing the rise of the defect density near the. surf ace. The 
Sah-Wang experiments [23] concerning defects from silicon-device 
processing suggest that vacancy complexes could give rise to the 
defects- observed, which is consistent with this value of r. 

Several considerations combine, however, to suggest that smaller 
values of r, and thus less complex configurations of the defects, 
will suffice to bring theoretical and empirical into conformity. 

First, if Auger recombination via.: traps preponderates in determining 


46 



18 3 

the recombination for a doping density of 2(10 )/cm-, then the 

needed value shifts to r<2, which suggests di-vacancies as a 
possible defect .configuration. Lack of experimental data about 
the pertinent -rate constants rules out a quantitative assessment 
of this possibility; the data available, however, indicate that 
band- to- band Auger processes play a second-order role in the 
device under 'study. Second, if the’ processing used to form 'the 
junction- were to depress x„(x) below its pre-processing value,, this 
also would drive down the needed value of the parameter r.. Finally, 
areal ihhomogeneity, discussed in Section III, would produce, the 
same tendency. 

Note that band-gap .shrinkage alone cannot explain the data. 
Even if the upper-bound estimate of its influence is assumed, a 
factor of 466 still separates experiment and theory. From a 
physical point of view, this occurs because for the dark cell 
under study essentially no holes are present in regions where the 
gap shrinkage might be appreciable. 

For this cell, subject to the validity of the assumptions 
made, the additional physical mechanism determining in this 

one-dimensional model for high-doping limitations is dominantly 
the drop in lifetime near the surface. 
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